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Abstract 

The deformation potentials of cubic semiconductors are re-examined from the point of view of 
the extended-basis sp 3 d 5 s* tight-binding model. Previous parametrizations had failed to account 
properly for trigonal deformations, even leading to incorrect sign of the acoustic component of the 
shear deformation potential d. The strain- induced shifts and splittings of the on-site energies of the 
p- and d-orbitals are shown to play a prominent role in obtaining satisfactory values of deformation 
potentials both at the zone center and zone extrema. The present approach results in excellent 
agreement with available experimental data and ab-initio calculations. 

PACS numbers: 71.15.Ap, 71.55.Eq, 73.21.La, 85.75.-d 
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The effect of uniaxial stress on the band structure of semiconductors has been a major 
theoretical and experimental topic for many years. With the development of strained-layer 
epitaxy, it has also become an important issue in modern material science and device physics. 
In their seminal approach, nearly half a century ago, Bir and Pikus established the strain 
Hamiltonian using the theory of invariants [1]. It depends on a number of deformation 
potentials describing the shifts and splittings of the various band extrema. For instance, for 
a given band near the Brillouin zone center, it reads as : 

W e = -a l (e xx + e yy + e zz ) -3b' 

Where are the components of the strain tensor, L is the angular momentum operator, 
{L x L y } = \{L y L x + L x L y ), and c.p refers to circular permutations with respect to the 
axes x, y, z. The coefficient a 1 is the hydrostatic deformation potential for the ith band, 
while b % and d % are respectively the tetragonal and rhombohedral (or trigonal) deformation 
potentials. We now explicit Eq. (1) for the T 6 and T 8 states of zinc-blende crystals : the T 6 
conduction-band energy only depends on the hydrostatic term owing of the L = matrix 
representation of the momentum operator. For the Tg valence-band edge, the heavy- and 
light-hole degeneracy is lifted (as L — 1) and the splitting depends on the strain orientation. 
Under [001] unixial stress or for lattice mismatched epilayers grown along the [001] direction, 
strain components can be written as: e xx = e yy ^ e zz and e xy = e yz = e zx = 0. Thus, the 
heavy- and light-hole bands split by an amount proportional to b. Under [111] stress or 
for mismatched epilayers grown along the [111] direction, we have e xx = e yy = e zz ^ and 
£xy = tyz = ^zxi^ 0, giving a valence-band splitting proportionnal to d. The two situations 
also differ by the presence or not of a static displacement of the anion and cation sublattices: 
in the former case all the atomic bonds in the strained crystal retain the same length and 
same angle with respect to the strain symmetry axis, giving no short range contribution 
to the strain hamiltonian. Conversely, for trigonal distortions, the equilibrium positions of 
atoms are no longer fully determined by stress invariants and a relative displacement of 
the sublattices is allowed. The resulting internal strain is represented by the Kleinmann 
parameter ( which ranges between and 1 [2|. The value ( = 1 corresponds a deformation 
with the same symmetry as the [111] strain but maintaining equal bond lengths of aoy/S/A, 
whereas ( = is related to the macroscopic strain that does not account for sublattice 
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d l [{L x L y } e xy + c.p] (1) 
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displacement. This gives simultaneously a long range or acoustic (tf) and a short range or 
optical (dg) contribution to the rhombohedral deformation potential d jjj: 



d = d'-jd (2) 

The strain hamiltonian has been extensively used in connection with the k.p theory to 
interpret experimental data and measure parameters. On a more fundamendal side, several 
ab initio calculations of the deformation potentials have been reported. More recently, atom- 
istic approaches using empirical parameters have become a strong challenger to k.p theory 
for a precise modeling of semiconductor nanostructures where compositions and deforma- 
tions can vary rapidly at the bond-length scale. A remarkable feature of atomistic theories 
(as opposed to the fundamentally perturbative character of the k.p theory) is their natural 
ability to treat the whole Brillouin zone. It follows that a "good" atomistic model must give 
a proper account of general distortions not only in the vicinity of the fundamental gap, but 
also at the edges of the Brillouin zone: the effects of strain actually is a very stringent test 
of atomistic models. Here we examine the effects of tetragonal and trigonal deformations 
from the point of view of the empirical tight binding (TB) theory. 

Within the tight-binding formalism, strain effects are mainly determined by scaling the 
Slater-Koster two-center integrals [4] (or transfer integrals) with respect to bond-length 
alterations, while bond-angle distortions are automatically incorporated via the phase fac- 
tors in the Slater-Koster matrix elements. This leaves a more than sufficient number of 
strain-dependent parameters to fit the deformation potentials at the Brillouin zone center. 
However, when trying to fit simultaneously the splitting of zone-edge conduction valleys, one 
encounters a difficulty. For the case of [001] uniaxial strain, it was shown in Ref. [6j that 
adding a term corresponding to a strain-induced splitting of the d-orbitals on-site energies 
(one-center integrals) leads to a much better overall fit of the deformation potentials at V 
and X. Indeed, since the Wannier functions of tight-binding models are Slater-type orbitals 
j^], the one-center integrals are expected to be sensitive to the environment of neighboring 
atoms. In principle, one should also introduce a [001] shear parameter of the on-site p en- 
ergies, and this can be generalized to all diagonal matrix elements in proportion to cubic 
and uniaxial distortions. However, in the atomic limit, the on-site properties should depend 
neither on strain, nor on chemical environment. This is nearly the case for the s and p 
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valence states that display an excellent degree of transferability. Conversely, it is clear that 
the excited s and <i-like orbitals have a strong free-electron character and corresponding 
energies must depend on strain-induced effects. Finally, one should remember that the "rule 
of the game" of atomistic models is to limit the number of empirical parameters to the min- 
imum required to account for symmetries and reproduce experimental (or ab initio ) band 
parameters within a given accuracy. For instance, the hydrostatic shift of on-site energies 
appeared to be useless parameters at the present level of model sophistication, as they are 
for a large part renormalized in the variation of transfer integrals with bond-length changes. 

To the best of our knowledge, the case of trigonal deformations has never been discussed 
in the framework of an advanced tight binding model. A first difficulty is the choice of a value 
for the internal-strain parameter (: contrarily to ab initio calculations, the atomic positions 
are an input of the tight binding model, not a result of the calculation. ( can be obtained 
theoretically either in ab initio calculations or from the fit of phonon dispersions, as first 
demonstrated by Nielsen and Martin 7 1. ( was precisely measured by X-ray diffraction forSi 
and Ge. The common value ( = 0.54 [8| is in agreement with first-principle calculations 
However, for GaAs, the most-cited experimental result £ = 0.76 [12] is still controversial 
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Inland dfifers significantly from the theoretical value C = 0.48 obtained by Nielsen and 
Martin [7J]. Note that the latter value gives an excellent representation of elastic constants 
and phonon frenquencies of GaAs and is corroborated by recent x-ray measurements which 
give £=0.55 [9]. The second difficulty is methodological since a fit of d comes out of a 
calculation, while there are two quantities to determine, d' and do. Here, d' is obtained by 
running the code using the fit parameters and setting ( =0, whereas a relative displacement 
of the anion and cation sublattices in absence of macroscopic strain, is used to calculate d$. 

The need for introducing a new shear parameter is evidenced by the failure of simpler 
tight-binding models : in the minimal sp 3 basis, the TB Hamiltonian cannot describe both 
b and it d satisfactorily, as seen for the diamond structure where b, d' , and do are linked by 
the analytical relations {3]]: 

16 

d = 16d' = T b (3) 
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For Ge, fitting &=-1.88 eV ±0.12 10] and considering ( = 0.54, one obtain: g?=-1.26 eV, in 
poor agreement with the experimental result: d=-5.0 ±0.5 eV [111]. The failure is mainly 
caused by an erroneous positive deformation potential d' (see Eq. 3) in sharp contrast with 



;he sign calculated by the self-consistent LMTO and ab initio pseudopotential approaches 
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14| which show a strong redistribution of the valence-electron density induced by the 



acoustic deformation. Similar discrepancies are obtained numerically for the zinc-blende 
semiconductors and no significant changes appear within the sp 3 s*d 5 approach, until a split- 
ting of on-site energies is introduced. The corresponding hamiltonian obviously depends on 
the strain direction. For a uniaxial stress along [001], the perturbation has the sym- 
metry, and the crystal field splits the fivefold degenerate d orbitals into two doublets and a 
singlet as : 



Edx Z = E dyz = E d (1 - <5 00 i(e 2Z - e xx )) 

Ed 2 2 = Ed 2 2 = E d (4) 

E dxy = E d (1 + 25 00 i(e zz — e xx )) 

For a uniaxial strain along [111] the perturbation has the symmetry and also splits 
the five equivalent d bands into two doublets and a singlet state. To handle this case, it is 
more convenient to rotate the coordinate system and choose the quantization axis z along 
the [111] direction. To avoid confusions we name "(111) basis" the new ( x, y, z) basis. 

The on-site d energies now corresponds to the representations Ai with E di _ 2 _ r2 , E x with 
Ed s - Z and E d , and E 2 with E d and E d : 



E d3i2 _ r2 =E d (l + 26 ul {e- z -e x )) 
E d „ = E d ,, = E d (l- 5 ul {e- z - e x )) (5) 
E d r, = Ed-- — Ed 

a ^2_y2 a xy a 

Negative e xy corresponds to conventional compressive stress. <5 oi and 5 m are shear 
parameters fitted to reproduce the tetragonal and trigonal deformation of the valence-band 
edge, respectively. Note that, since perturbations result from different modifications of the 
nearest neighbors positions, there is no reason for an exact geometrical relation linking 5ooi 
and Sin. The scheme of level splittings for uniaxial compressions along [001] and [111] 
is shown in figure 1. Following this procedure, we demonstrate excellent agreement with 
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experiment for b and d. In addition, the figures coming out for the acoustic deformation 

n 

potential d' are consistent with self-consistent LMTO results Il3fl . For instance, for Ge, 

n c 

using the parametrization of Ref. [6J, with £ oi=0.54 and 5 m =-1.5 , we get: 6 =-1.9 eV, d 



lOj, d=-5.0 



=-4.6 eV, and d' =-1.37 eV, in agreement with experiment: &=-1.88 ±0.12 eV 
eV ±0.5 Q, and the LMTO calculation: d'= -1.3 eV Q- 

Furthermore, as discussed in Ref. f| , introduction of <5ooi allows to obtain simultaneously 
a fit of b and of the conduction-band splitting of X valleys under [001] stress. However, for 
[111] stress, when examining the shear deformation potential of the L\ conduction extrema, 
using our fit value of 5m, we find a very large and disapointing discrepancy: D\ c =12.h 
eV, to be compared to the experimental value of 18.3 eV llj. This witnesses that there is 
still a missing parameter to describe properly the effects of trigonal deformations. From a 
quantum chemistry point of view, this discrepancy can be understood from the fact that the 
L-conduction band minimum is dominated by s and p states, in contrast with wavefunctions 
at X that have a strong d character [6j. Therefore, the splitting of on-site d energies 
cannot help much to fit the strain deformation potential at L. Instead, the splitting of 
on-site p energies, which was up to now a non-necessary parameter, becomes important to 
model strain field anisotropies. At this point, the internal logic of the model implies that if 
the strain induced splittings of the one-center integrals play such an important role, their 
shift under hydrostatic strain should also be taken explicitely into account instead of being 
renormalized in a strain dependency of two-center integrals. Using again the "(111) basis", 
the corresponding contributions to (lll)-strain hamiltonian are written as: 



e p* = E Py = E p (! - - e 2 )) 

E p ^ = E p (l + 27T lu (e- z -e s )) 



(6) 



The values of d, d', do, and D\ c , for Ge and GaAs semiconductors are compared in 



table 1 with available experimental data and LMTO [13j evaluations. We use values of 
( issued from experiments and fit 7r m and 5m to reproduce the trigonal deformation of 
the valence-band edge and L-conduction band valleys. To achieve a complete description 
of strain effects in the spds* model, on-site Hamiltonian matrix elements were also scaled 
with respect to bond-length changes resulting. Althgouh the number of tight-binding 
parameters is increased with respect of Ref [6j their numerical determination through 



a multi-parameter fitting procedure still converges very well and resulting values gain 
a more intuitive view of chemical dependences. A detailed description of the complete 
parametrization is beyond the scope of the present work and will be the subject of a 
forthcoming paper. The data presented in Table 1 demonstrate that our results are in good 
agreement with experimental and theorerical values for d an D\ c . The magnitude and sign 
of d' and do are well reproduced compared to the LMTO calculations. Our values for do 
differ from the experimental results whoses values are actually contreversial, a discrepancy 
observed already earlier; (compare Table 1). 

In conclusion, we have demonstrated that the sp 3 d 5 s* model requires diagonal matrix 
element shifts to correctly reproduce uniaxial [klm] strain for cubic semiconductors. In 
order to test our model we have calculated the acoustic and optical contributions to the 
trigonal deformation potentials and found a good agreement with experiment and LMTO 
results. A major improvement compared to smaller TB models was the correct sign and 
magnitude of the acoustic deformation potential d' directly related to the shear parameter 
of d-states. This TB model provides a valid framework for the calculation of strain effects 
in self-assembled quantum dots. 
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FIG. 1: Schematic plot of energy-level splitting of d states induced by an uniaxial strain along 
the [111] (left panel) and [001] directions (right panel). For [111] strain, the quantization axis z is 
chosen along the [111] direction. 

TABLE I: Comparison of trigonal deformation potentials obtained in present work with LMTO 
and experimental results. We use £ = 0.54 for Ge and Q = 0.55 for GaAs. 



Ge GaAs 
d d! d D\ c d d' d D\ c 
this work 30 -0.9 -5.0 17.1 27 -0.8 -4.5 18.0 
LMTO a 22.4 -1.3 -5.0 25 -0.99 
Expt. 6 M -5.0 18.4 44 -4.5_20.8 



a Ref. [13|], b Landolt-B ornstein Ref. [ill]. 
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